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ABSTRACT 

Stabilizing a nulling interferometer at a nanometric level is the key issue to obtain deep null depths. The 
PERSEE breadboard has been designed to study and optimize the operation of a cophased nulling bench in the 
^> most realistic disturbing environment of a space mission. This presentation focuses on the current results of the 
HH PERSEE bench. In terms of metrology, we cophased at 0.33 nm rms for the piston and 80 mas rms for the 
f~| tip/tilt (0.14% of the Airy disk). A Linear Quadratic Gaussian (LQG) control coupled with an unsupervised 
Oh vibration identification allows us to maintain that level of correction, even with characteristic vibrations of 
Q , nulling interferometry space missions. These performances, with an accurate design and alignment of the bench, 
i—^ • currently lead to a polychromatic unpolarised null depth of 8.9x 10~^ stabilized at 3x 10~^ on the [1.65 — 2.45] |^m 
c/2 , spectral band (37% bandwidth). 
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> ■ 1. INTRODUCTION 

ff^ IR spectrometry is a powerful method to characterize the hundreds of exoplanets discovered by indirect methods, 
■ and ultimately to search for life on telluric ones, as it allows to evidence atmospheric signatures of interest (H2O, 
, CO2, CH4,...) boosted by the thermal emission of the planet. The main challenge is to isolate the planet 
0^ ' from its parent nearby star; one of the very few possible methods is spaceborne nulling interferometry, allowing 
. implementing very efficient coronagraphs at the diffraction limit. Within this framework, several space projects 
such as DARWII^ilE TPF-PS, FKSP or PegastP^I^ have been proposed in the past years. Most of them are 
based on a free-fiying concept to create hectometric pupils with an acceptable launch mass. But the price to 
pay is to very precisely position and stabilize the spacecrafts relatively to each other and also with respect to 
^> [ an inertial direction, in order to achieve a deep and stable extinction of the star. Understanding and mastering 
5^ ' all these requirements is a difficult challenge and a key issue towards the feasibility of these missions. Thus, we 
decided to experimentally study this question and focus on some possible simplifications of the concept. 

Since 2006, PERSEE (Pegase Experiment for Research and Stabilization of Extreme Extinction) laboratory 
test bench is under development by a consortium composed of Centre National d'Etudes Spatiales (CNES), 
Institut d'Astrophysique Spatiale (IAS), Observatoire de Paris-Meudon (LESIA), Observatoire de la Cote d'Azur 
(OCA), Office National d'Etudes et de Recherches Aerospatiales (Onera), and Thales Alenia Space (TAS)®. It 
is mainly funded by CNES R&D. PERSEE couples an infrared wide band nulling interferometer, based on a 
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Modified Mach-Zehneder (MMZ) and a geometrical achromatic phase shifter, a disturbance simulator and a 
cophasing system based on a LQG tracker for the optical path difference (OPD). The first results of the bench 
were already presented in Ref. [HI in this paper, we focus on complete it with the latest improvements and results 
of the bench. 

After a short description of PERSEE goals and setup (focused on recent changes), we will present in Sec. |3] 
its various calibration procedures. Recent results are reported in Sec. 2] for the cophasing system and in Sec. [S] 
for the nuUer. 

2. PERSEE DESCRIPTION AND UPDATES 

This section provides the reader with a very brief description of PERSEE. Much more details can be obtained 
by the interested reader from Refs. [51 ITTHH1 

2.1 PERSEE goals 

The goal of PERSEE is not to reach the deepest possible nulling. Starting from a state of the art nuUer of the 
2006-2007 period, it is an experimental attempt to better master the system flowdown of the nulling requirements 
both at payload (instrument main optical bench) and platform levels (satellites) . The balance of the constraints 
between those two levels is a key issue. The more disturbances the payload can face, the simpler the platforms 
are, the lower the cost. The general idea is hence to simplify as much as possible the global design and reduce 
the costs of a possible future space mission. 

The detailed objectives have been described in Ref. 8. Our main requirement is to reach a 10~^ nulling ratio 
in the [1.65-2.45] |^m band (37% spectral bandwidth) with a 10~'^ stability over a 10 h time scale. Another 
important requirement is to be able to find and stabilize fringes which have an initial drift speed (as seen from 
the interferometer core) up to 150 |j.m.s"^, as this can greatly simplify the relative metrology and control needs. 
We want to study and maximize the rejection of external disturbances introduced at relevant degrees of freedom 
of the optical setup by a disturbance module simulating various environments coming from the platform level. 

2.2 Modifications of the source module 

In the original design, three sources were used: a laser diode at 830 nm and a SLED at 1320 nm for the dual band 
metrology system, and a black body used in the band [1.65 — 2.45] |J.m for the science path. The three sources 
were coupled in single mode optical fibers, that were joined close to each other at the focus of the collimator, 
to simulate the star. That configuration was not efficient, because the SNR in the camera was too low, the 
science spectrum was polluted by the SLED fiux, and the contrast was too low in the Fringe Sensor (FS). So we 
had to change the configuration at the focus of the collimator, to have only one fiber that simulates the star. 
Unfortunately, the black body was not powerful enough for the FS detectors to be used for metrology in addition 
to science. 

The adopted solution was to use an infrared supercontinuum source, a polychromatic source powerful enough 
to emit both in metrology and science paths. Its pulse rate is at 20 MHz, much higher than the loop frequency 
of the FS (1 kHz) and the camera (100 Hz). It was not possible to connectorize the supercontinuum fiber, so 
we designed a simple injection bench with two off axis parabola, to inject the fiux in a single mode PCF fiber. 
The flux was even too high for the science camera, so we have to defocalize the PCF fiber to reduce the coupled 
light. 

Another modification of the source module, presented in that conference,'^^ will include a planet and a 
exozodical disk. 



2.3 Flux reference for the nulling rate analysis 

On the science spectrometer, two flux are measured in order to perform tfie calculation of the nulhng rate on 
each spectral channel i, by 

^ (1) 

where Imin is the flux of the destructive interference, directly obtained at the nulled output III of the 
MMZ. Imax is the flux of the same output, in the case of a constructive interference. It cannot be measured 
simultaneously with Imin , so we firstly used the flux of the related constructive interferometric output II of the 
MMZ. This solution has two disadvantages: firstly, the single-mode fibers of outputs II and III must be perfectly 
aligned, in order to have optimal injection in both outputs, and the two fibers should not be subject to differential 
displacements. Secondly, the output II is a coherent fiux, so it is sensitive to OPD variations. Especially in that 
case, the differential phase at the output II is not zero but 21 ° , thus the output II is not strictly constructive 
and does not measure exactly Imax- 

For those reasons, we modified that reference channel, by positioning the reference fiber in front of the Mq 
collimator, selecting a part of the unused flux of the source, before the interferometric train. 

The flux Iref of that output is not equal to the flux Imax, but they are proportional: Iref = jimax- Also, 
the dynamics of the camera (« 10^) is not sufficient for the direct measurement of the desired nulling rate, that 
is why we use an optical density of transmission T when we calibrate the coefficient 7, and we remove it for the 
measurement of Imin ■ The goal is to reduce the contrast of the measurement to fit in the range of the camera. 
Finally, Eq. ([1]) becomes 

AT — Imin,i 

' ~ T I t ' ^ ' 

3. CALIBRATION SCHEME OF THE BENCH 

PERSEE bench has two control loops in piston and tip/tilt, which require a few calibrations, as well as the 
scientific path. There are three types of calibrations. The first type (Sec. 13. 1113. 3p depends only on the configu- 
ration, and need to be applied just once. The second type of calibrations (Sec. 13.4113.61) depends slightly on the 
environment, especially thermal drift. The last type of calibrations fSec. 13.71 13. 8p is strongly affected by thermal 
drift. 

3.1 Calibration of the spectral bands 

To calibrate the spectral bands, we used the method of Fourier Transform Spectroscopy (FTS). The FTS uses 
the properties of the Fourier transform to determine the spectrum of a source with an interferometer. In most 
cases, the interferometer is a Michelson with a mirror mounted on a translation, to scan the interference fringes. 
With a MMZ, the situation is slightly different. For each wave number a, the intensity as a function of S is 

1(6, cr) = /(cr)[l + V{a) COs(27rcr(5 + (/)(cr))], (3) 



with I — la + lb and V — 2^ I ah '/ {la + lb)- I (c) is the spectrum that we want to determine. Then the 
intensity measured at the output of the interferometer is 

POO 

I{S)= / I{S,a)d<j. (4) 
Then the Fourier transform of I (S), I (cr), is 



/(a)=Vo(a) + -e^^('^)/(a)F(a), 



(5) 




Figure 1. Left: Spectrum of the two spectral bands of the FS, which have two spectral bands over 4 outputs. Right: 
Spectra of the nine channels of the science camera, each one spread on two pixels on the detector. 



with So the Dirac delta function, and /o = /q°° /(<T)d(T the total flux. We find the case of the Michelson 
interferometer, when the contrast is maximum and the phase is zero: I{<t) = IqSo{<j) + ^/(cr). In the case of the 
MMZ, we do not measure the spectrum of the source, but the product of the spectrum and the contrast. 

Figure [T] presents the spectra of the FS bands and the science path. On the FS spectrum, we can see the 
pump laser of the supercontinuum source, at 1064 nm, and the absorption of the water vapor around 1.4 |j.m. 
The four outputs of each FS spectral band are slightly different. That effect deteriorate the linearity of the FS, 
especially on the J band. On the camera spectrum, we can see that the different channels are not equal, in 
intensity and bandwidth. The absorption of water vapor around 1.8 |^m is also visible. With this analysis, we 
can deduce the central wavelength and the bandwidth of each channel. 

3.2 Calibration of the density 

We have three optical densities on the scientific spectral band, but they were not calibrated by the manufacturer, 
who did not have the appropriate instruments in this spectral band. The theoretical absorbance of those densities 
are 3, 4 and 5. So the idea is to calibrate the densities together, to reduce the contrast between measurements. 
So we compared density 3 with 4, 3 with 5, and 3+4 with 5. The needed contrast is then 10^, which is much 
less than the camera limit of 10''. 

3.3 Calibration of the actuators 

It consists in calibrating the PZT actuators of the two active Mg mirrors. Each one has three piezoelectric 
motors, which allow movements in piston and tip/tilt. This calibration process measures the interaction matrix 
and computes its generalized inverse to derive the command matrix. The interaction matrix Mint is obtained 
by sending known commands (voltage ramps) to each piezo motor, and by measuring the slope on each axis on 
the Field Relative Angle Sensor (FRAS). The matrix Mint is then defined by 




M,;. 



dV2 



(6) 



We perform then a Singular Value Decomposition (SVD) to derive the inverse. The piston mode is not 
seen by the FRAS but can be reconstructed from a geometric model of the actuators, and its control added in 
the command matrix, which transforms commands in piston and tip/tilt to voltages to be sent on the 3 piezos 
120 ° apart in azimuth. 
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Figure 2. Reconstruction of the convolution of the PSF and the fundamental guidance mode of the single-mode fiber, for 
each spectral channel of the science camera. 



3.4 Calibration of the tip/tilt reference position 

The purpose of this calibration is to obtain the reference position of the Mg mirrors to inject a maximum flux in 
the single-mode fiber that guides the flux from the output III of the MMZ to the scientific camera. A sequential 
scan of each mirror is performed to move the focal spot in front of the fiber; the signal is synchronously acquired 
and is used to reconstruct an image. The theoretical intensity can be approximated by a gaussian functiocP^. 



tnj,max 



exp 
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(7) 



with linj.max the maximum injected flux, D the diameter of the beam, a and /3 the angles of tip and tilt, 
and aref and /3re/ the reference position which maximizes the injected flux. 

The result of that reconstruction is shown in Fig. [2l for each channel of the science camera. We get a spot 
corresponding to the convolution of the PSF for each spectral channel with the fundamental guidance mode of 
the single-mode fiber. According to Eq. ([7]), the diameter of the spot is proportionnal to A. This can be seen in 
Fig.O the diameter of the spot decreases as the wavelength decreases (left to right). We compute the centers of 
gravity for the two arms, which give the offsets to be applied to both mirrors. As the spots are centro-symmetric, 
this position corresponds to the maximum injection flow of each arm in the fiber. The dispersion of the centers 
of gravity for the different channels are negligible, so there is no chromatism in the reference position in tip/tilt. 

3.5 Calibration of the fringe sensor 

That procedure has already been described in Ref. 1131 Since then it encountered a few changes, but the idea 
remains the same: we inject an OPD ramp with the Mg actuators, to measure the phase shifts between the 
4 outputs of the MMZ on the two spectral bands. Then we measure sequentially the flux on each arm, and the 
dark current of each detector. We have then the matrix Mj^s given by the relation 
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with la and lb are the flux of the two arms and /iEnv((5) the fringes envelope. The matrix Mps is inverted, 
to obtain the demodulation matrix D. 
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Figure 3. Interferogram acquired during the calibration of the camera. 



3.6 Calibration of the science camera 

The calibration of the camera is made along with the calibration of the FS. Therefore, it uses the same sequence: 
an OPD ramp, the flux in each arm, and the dark current. With those information, we can analyze a few 
parameters, and compute the coefficient 7 of Eq. ([2]) for the calculation of the null rate. 

With the data flux of arms a and b, la and lb, we can estimate the photometric unbalance e by 



{h), + {ia); 



(9) 



As 

Imax — (1 "h V^)(/a ~h Ih) i with. the COntiaSt \^ — \^ then Imax 

~ 2(/a + /fe); we estimate the coefficient 7i 

for each spectral channel by 

Figure 13] presents an interferogram obtained on the camera during the OPD ramp of the calibration. We can 
see that every channel has different wavelengths, but their minima are at the same position. With those data, 
we can use a FTS mode, to determine the central frequency fi and the phase 0^ of each sinusoid. 

With the frequency fi, we can determine the wave number ai with the slope of the OPD ramp v: (Ti = fi/v. 
With the phase (pi, we can determine the position of the minimum of each sinusoid, by 

(11) 



— Ore/ -1- 7; , 

iTTCFi 

with 5ref is the reference position of the FS. 
3.7 Correction of the tip/tilt reference position 

When the temperature drifts, the reference position of the FRAS also drifts. This shift is usually small, so it is 
not necessary to repeat the procedure described above. A more efficient procedure has been developed for this 
purpose. It is inspired by the procedure of correction of the OPD position described below, called ditherin^iSl_ 

Close to the reference position, we can approximate the exponential term of Eq. ([7|) with a quadratic term: 



/(a, P) ~ hnj,max {l " y [(" " + (/? " Preff] ) (12) 

For the two arms, for both axis, we put a smah angle 9 — Q arcsec (11% of the Any disk), and we measure 
the flux injected in the fiber. If ao and /?o are the current positions, we note /q = ^(o^Oi /^o)i ~ ^(o^o ~ 
Ia+ = /(ao + 0,i3o), = I{aQ,Po - 9) and Ifi^ = /(ao,/3o + 6). 

Then, we calculate the reference position with 



aref = ao 



Pref = /3r ^ 



2 + /a_ - 2/o) 



(13) 



2 (//3+ + Ip_ - 2/o) 



This procedure was successfully implemented on the bench, and showed a very good accuracy (« 50 mas, 
0.1% of the Airy disk). 

3.8 Correction of the OPD reference position 

Despite the calibration procedure of the camera described above, the accuracy is not enough to guaranty a sub- 
nanometer positioning of the OPD to minimize the null rate. That position is also very sensitive to temperature 
variations. To correct the reference position of the FS Sref, we use the technique of ditherin^i^. 

The null rate has a quadratic dependence with the OPE^l 

N{d)^N„,,n + {7:a{d-6ref))^ (14) 

So we measure the null rate for three positions of OPD: the current position Sq, and two other positions Sq — e 
and So + e, with e = 2 nm. If we note A^o = N{6o), N^^ — N{So~e) and A^+e = N{do + s), the reference position 
is then 



2 (A^+, + N^, ~ 2No) 



That reference position is very sensitive to temperature variations, about 600 nm.K"^. However, we consider 
that temperature variations over 100 s are less than 1 mK in normal conditions, so the thermal drift is less than 
0.6 nm over 100 s. 



4. LQG EXPERIMENTAL VALIDATION 

PERSEE aims to simulate a nulling space mission, constituted of three satellites: two siderostats and a central 
combiner. The three satellites are controlled by impulsive gas thrusters and reaction wheels, from individual 
star trackers, a radiofrequency metrology, optical lateral metrology and off-loading signals from the central 
combiner. Those systems lead to different types of disturbances, viz., microvibrations due to the reaction wheels, 
linear-parabolic drift due to the solar pressure and the gas impulses, and low-frequency pointing control (PC) 
residue. 

Vibrations are the main disrupters of the OPD, which are not always well corrected by an integrator in a 
control loop. We implemented the method described in Ref. 13, based on the work for the extreme adaptive 
optics instrument SPHERE at Onera^^. It is based on a Linear Quadratic Gaussian (LQG) control law that 
provides an optimal correction of vibrations, in the sense of residual phase minimum variance. 



POINTING CONTROL RESIDUE CORRECTION 
Controller c c!c,pc PcPC 
Integrator 0.45 nm rms 0.17% 
LQG 0.37 nm rms 0.09% 
Table 1. Correction of the PC residue with the integrator and the LQG. 



4.1 Identification of the pointing control residue and vibration filtering efficiency 

In a first test, we analyze the correction of the PC residue, with a classical integrator and a LQG. 

We note CTq^ ^ the open- loop OPD variance for a case where we inject a PC residue of variance af^^ , and 
vibrations of variance crf„^^j^. If (JoL,exp is the open-loop OPD variance of the experiment without any injected 
disturbance, then 

'^OL.dis — ^OL.exp + '^inj,PC' + '^inj.vib- (16) 



So if ctq^, PC is the open-loop OPD variance for the case where we inject only PC residue, we have CTq^, pc = 

'.,exp + '^inj,PC^ ^'^'i i'^^'^ '^OLMis — ^OL,PC + ^'inj,vib- 

For this test, we compute the following power residue percentage 



^2 ^2 

Pc,pc = 2^^'^ _ "'l'^'^ , c e {int,LQG}, (17) 

'^OL.PC '^c,exp 

with cr^ the residual OPD variance for the controller c (integrator or LQG) when we do not inject distur- 
bances, and cr^ pQ the same quantity for the case where we inject the PC residue. In that case, if ct^ — a'^ ,,^p, 
then Pc^pc = 0%, the injected disturbance is totaly corrected, and if <J-^ pc = <Jq]^ p^, then Pc,pc = 100%, the 
controller has no effect. 

Numerical results of the two controllers are presented in Tab. [TJ With that type of disturbance, i.e. with a 
continuously decreasing spectrum like in AO with turbulence, the integrator and the LQG have similar perfor- 
mance: the power residue is very low in each case, with a slight improvement from the LQG. 

For the following tests, we will compute a vibration residue percentage, similar to the one used in Ref. 1191 
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ce{int,LQG} (18) 



with tTcdis is the residual OPD variance for the controller c for the case where we inject the PC residue 
with vibrations. So if cr'^j^^g = cr^pc^ the pi — 0%, the vibrations are totally removed. On the contrary, if 



c,dis 

(jfni vib' Pc — 100%: the vibrations are not corrected. 



4.2 Experimental correction of a single vibration 

In those tests, we inject the PC residue with one vibration, at several frequencies, and we compare the correction 
of the integrator and the LQG. 

Table [2] presents the results of the integrator and the LQG for different vibration frequencies. As expected, 
the correction of the integrator deteriorates as the frequency increases, in contrast to LQG, which maintains a 
fairly constant performance plqg close to 0%. 

For fyib ~ 1.1 and 2.2 Hz, the injected vibration was not identified, but corrected as a part of the PC residue 
contribution (the continuum part); it does not affect the result of the correction: the disturbance is totally 
removed. 



CORRECTION OF A SINGLE VIBRATION 
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Figure 4. Frequency evolution of the identification, and on the residual OPD (simulation and experimental results). 



Table [2] shows that after 22 Hz, for 10 nm vibrations, the integrator is not efficient to correct the disturbances 
at the required level of 1.8 nm rms. Also, the power of the vibration at 110 Hz is doubled, because of the 
overshoot in the integrator closed-loop transfer function. 

These tests revealed a limitation of identification: the spectral resolution of the FFT. Indeed, if a sinusoid 
has not an integer number of periods in the signal supplied to the identification, the fitted vibration will have a 
damping factor greater than reality. 

Figure m illustrates that effect. In that example, the model used is a typical PC residue, with a single vibration 
around 22 Hz (black solid line). For the estimation 1 (blue solid line), the vibration has an integer number of 
periods in the time sequence used for identification. The estimation is then close to the model for the PC residue, 
and the identified damping coefficient for the vibration is 10~^. Then the OPD residue is close to the ideal one, 
giving almost a white noise. For the estimation 2 (red solid line), the number of periods is not aninteger, the 
peak maximum is between two points of the periodogram. In that case, because of the contrast between the 
vibration and the PC residue, the peak is not sharp, but has some leakages around the central frequency. The 
identification is degraded, the PC residue is not perfectly fitted, especially at the cutoff frequency. For the 
vibration, the identification gives a greater damping coefficient of 5 x 10"'*, and a reduced amplitude. On the 
OPD residue, we can see that the signal is overcorrected around the vibration frequency, while the vibration is 
not perfectly suppressed. This overcorrection is compensated by a higher residue at low frequencies. 

4.3 Experimental correction of realistic disturbances 

These tests aim to simulate as accurately as possible the conditions of disturbance of a space mission. To do so, 
we consider two levels of PC residue: an optimistic value of 3 nm rms, and a conservative value of 15 nm rms. 




Figure 5. PSD and cumulative OPD for the correction of realistic disturbance, with an integrator and a LQG. The 
parameters of the disturbance are frm = 4.5 Hz and (tol,pc ~ 15 nm rms. 

COMPARISON B/W SIMULATIONS AND EXPERIMENTS 
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Table 3. Influence of the number of vibrations. 



Added to that residue, we also inject the characteristic vibrations of 6 reaction wheels, rotating at different 
frequencies, close to a central frequency frw We consider three rotation frequencies, 1.5, 3 and 4.5 Hz, corre- 
sponding to typical GNC design values. Those reaction wheels generate a set of vibrations, related to the first 
harmonic of the wheel. 

No detailed mechanical model was available, and the following assumptions were made from typical satellite 
behavior known at CNES. Some of the frequencies are amplified by structure modes in the satellites: a mode of 
sunshade at 15 Hz, a fundamental lateral mode of the satellite at 55 Hz, and local modes of the wheel cavity 
and the optics at 75 and 85 Hz. 

Those parameters lead to six different cases that are studied with a simple simulator of the bench, and 
compared to experimental results. For a central frequency frw = 1-5 Hz, only the first harmonics are significant, 
no harmonic is amplified by the structure modes. For frw = 3 Hz, some harmonics are amplified at 55 Hz, while 
for frw = 4.5 Hz, harmonics are also amplified at 75 and 85 Hz. 

Figure O presents the correction of the worst studied case: frw — 4.5 Hz and (Jol,pc = 15 nm rms. For the 
LQG, the identification is made on the 20 strongest vibrations. We set the integrator gain at 0.5. 

Table [3] presents the comparison between the simulations and the experimental results for the 6 different 
cases. We can see that the integrator can manage the disturbance at the required levels for frw = 1-5 Hz, 
whatever the level of RF PC residue. For those cases, the LQG gives a substantial improvement, leading to a 
nulling contribution divided by 2. The interest of LQG is more significant for the other cases, because amplified 
harmonics are not properly corrected by the integrator. As expected, the LQG corrector removes almost the 
whole power of the disturbance. 



5. NULLING ANALYSIS 



5.1 Nulling calculation and contributors 

On PERSEE, we carry out a spectral analysis of the null rate, so we analyze separately the different spectral 
channels. Thus, we calculate an average null rate over the entire band, rather than an integrated null rate, which 
is more dependant of the spectrum. Therefore, the spectral bandwidth is larger. 

With the result of Fig. [1] we calculate the spectral widths for each channel i. Then the total spectral 

bandwidth is deduced from the extreme channels: 



Ao-„ , A(Ti 

Af7„i = CTn H CTi 



2 ' 



(19) 



with n the number of channel {n = 9). Then, the total bandwidth is — 37% around the central wave 

number cTm = 0.5 |^m"^, i.e. a central wavelength Am — 2.0 |xm. For the analysis of the bench performance, we 
will then compute the mean null rate — 

The null rate can be decomposed in a sum of independent contributorJi^ 



(20) 



where: 

• A<j)c{t) is the phase fluctuation at the central wavelength, while Aip\{t) takes into account the chromatic 
dispersion of the null; 

• Odia is the angular diameter of the source, amax the maximum wavenumber of the spectrum and b the 
interferometric baseline; 

• A(/)s-p is the phase difference between the two polarizations s et p, and a^ot is the rotation angle of the 
polarization s (or p) between the two beams of the interferometer; 

• e{t) are the fluctuations of intensity mismatch, defined by e = , with la and the flux on the arms 
a and b. 
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Figure 6. Best result of null rate on 100 s: temporal sequence of the mean null rate (left) and temporal average of the 
null rates on the different spectral channels (right). 



Figure [S] presents the best result obtained on the bench, on a 100 s timescale, with an acquisition frequency 
of 100 Hz, and an integration time of 2.58 ms. The mean value is 



(^m)t=ioo. = 8.9 X 10-6. 

The standard deviation aN,T depends on the integration time T. We chose to consider T = 1 s, more 
representative of a real mission. Then, we have crAr,T=2.58ms = 1-8 x lO^^ and <tn.t=1s — 2.7 x lO"''. The 
chromatic dispersion of the null rate is between 6.9 x 10"^ and 1.58 x 10~^ on the spectrum, the degradation of 
the null rate at short wavelengths is explained by both OPD and flux mismatch effects, but the degradation of 
the highest wavelength is probably due to polarization effects. 

5.2 OPD fluctuations and chromatism 

The OPD fluctuations are an important contribution of the null rate. The different spectral channels are disturbed 
differently by the OPD fluctuations. 

Figure [7] presents the contribution of the OPD on the mean null rate. Compared to Fig.[6l we can see exactly 
the same disturbances, with a similar amplitude. The OPD disturbances are characterized by their standard 
deviation aoPD — 0.80 nm rms. According to Eq. [201 the contribution of the OPD on the null rate is then 



Ns = ^2 {af)^ {{S - 5„0'), , with <5„, = (21) 

That contribution can be reduced by having Sref — Sm, especially with the calibration described in Sec. 13.81 
But for this case, the calibration was not operationnal, so the search of the optimal dref was made manually. 

So the contribution of OPD fluctuations and static error are 



Ns^7r^a^)^{Slf. + alpj,) (22) 

In this case, we measured \Sref \ — 0.48 nm. So the OPD contribution is Ns = 1.7x 10~®. We can also calculate 
the contribution of the null rate stability <jn,s.t- we have (TN^T=2.58ms = 1-8 x 10~^ and <jn,t=1s = 3.1 x 10"''. 
Those values are almost equal to the stability of the null rate, which shows that the null rate is only disturbed 
by the OPD variations. 

The reference position of the FS does not optimize the null rate of each spectral channel, because of the phase 
dispersion (A0>,(i))j. That dispersion is measured during the correction of the OPD reference position presented 
in Sec. [31 




Figure 8. Chromatic dispersion of the positions of best null rates for the channels of the camera, and its contribution to 
the null rate. 



The contribution of the chromatism to the null rate, is therP^ 



(23) 



So the total contribution on the whole spectral band is 



(24) 



Figure[5]presents a measurement of the chromatic dispersion on the channels of the camera. The contribution 
of the chromatism, calculated with Eq. ([M)) . is then N\ = 4.6 x 10~^. This chromatism is due to the beam 
splitters of the MMZ, especially the thickness differences of both plates and coatings, and also the difference in 
incidence angle between the two arms. 

5.3 Stellar leakage 

The contributor corresponding to stellar leakage N^, , corresponds in Eq. (PP)) to 



The angular diameter 9ciia of the source is given by the ratio of the diameter of the fiber d and the focal 
length /o of the collimator. On PERSEE, d = 4.3 [a.m and /o = 750 mm, 6 = 50 mm and dmax = 0.61 M-m"^, So 
the stellar leakage would be iV* — 1.9 x 10^^. 

However, we obtained much better null rate. Indeed, this calculation applies only for a spatially incoherent 
source, whereas in the case of PERSEE, the source is spatially coherent, due to the single-mode fiber. Therefore, 
there is no stellar leakage. 

5.4 Photometric mismatch 

The contribution of the photometric mismatch depends on the tip/tilt residue of the control- loop, and the static 
photometric mismatch of the instrument. 

From the Equation (IT^ . we can write 



N^: = —{Odia<^i 



max 



hf. 



(25) 




(26) 
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Figure 9. Photometric mismatch on the different channels of the camera. 

with Cl = {aa - 0.a,refY + {Pa " fia,reff and Q = {ab - Ub^reff + {Pb " Pb^ref)'^ ■ 

Then the photometric mismatch on each channel i is 

er = e,tat,^ + ^<y^D\C-eb), (27) 
with Sstat.i — ^^'^^"^^^ ^a,i,max ^ ^^iQ Contribution of photometric mismatch on each channel is 

Ne,= '^ + ^^<jtD\C~a?- (28) 
The contribution on the null rate is then 



The static contribution of the photometric mismatch £stat,i is calculated during the calibration of the camera, 
and is only due to the instrument optics. Figure IHl presents a measurement of Estat, and its contribution on the 
null rate, calculated by 

NeM = (30) 

Thus, the total contribution of instrumental photometric mismatch is given by 

Ne,stat = 1^^. (31) 

That contribution is then N^^stat — 3.9 x 10^^. 

The residue of tip/tilt control is also measured, and has a typical standard deviation Uup/ tut, axis = 56 mas rms 
on each axis, so the stability of the position on each arm is <Jtip/tut = 80 mas rms (0.14% of the Airy disk). 

According to Eq. (0^ . the dynamic contribution of the control loop residue to the null rate is 

Ne,dyn^^D^af)^at,^/,,^,. (32) 



Then, that contribution is N^^dyn — 1-1 x 10 ■ It is clearly negligible compared to Ng^stat, but the relation 
™ '^tip/tut very sensitive to the FRAS correction residue. 

For the stability over 100 s, the contribution of the tip/tilt control loop is very faint, at 'JN.tip/tut,T=is = 

1.4 X 10^^^, a few decades under the contribution of the OPD control loop. 

5.5 Polarization study 

We have no way to measure the polarization contributions, but we can estimate them from subtracting the other 
known contributions, by the equation: 

Npoi = N^~Ns-Nx~N, (33) 

So, with Nm = 8.9 x 10-6, Ns = 1.7 x 10"^, A^a = 4.6 x lO"'^ and = N^^stat = 3.9 x lO''^, we can estimate 
the contribution of polarization effects at Npoi — 2.8 x 10^^. 

So if that contribution is only attributed to the phase difference between the polarizations s and p, the OPD 
difference between the two polarizations is S^-p — 2.1 nm. On the contrary, if the contribution is only comming 
from the angle of the polarizations of the two arms, that angle is arot = 3.3 mrad. 

5.6 Stability of the null rate over 7 hours 

We are currently conducting stability tests over a few hours (typically 10 h), the typical duration of a target 
observation during a space mission. The first test was made over 7 h during the night, with a sub-optimal 
alignment. The average null rate is then higher than the value analyzed above. 

Figure fTOl presents the first long-term measurement of the null rate over 100 s. We can see on the left figure 
that the nuU rate has a higher value than the best result: {Nm)rp^^^ = 3.37 x 10~^, but shows a very good 
stability over 7 h: ctat.t^ioos = 9 x 10"''. That measured stability is a decade better than the specification: 

<JN,T=Was — 10^ . 

In Fig. [TO] (right), the null rate spectral profile looks like the one presented in Fig. [SJ the nuU rate is better 
at the center of the spectrum than at the edges. The contribution of OPD and flux mismatch are also higher at 
shorter wavelengths, and the contribution of polarization should be higher at longer wavelengths. The chromatic 
dispersion of the null rate is between 1.3 x 10^^ and 8.7 x 10^^, so the amplitude is 7 times higher than on the 
best result. 

Table m summarizes the null rate analysis developped in this section. In the near futur, we should be able to 
achieve an average null rate of 10~^ over 10 h, with a stabihty better than 10^^, i.e. a decade under the original 
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Figure 10. First measurement of null rate over 7 h: temporal sequence of the mean null rate (left) and temporal average 
of the null rates on the different spectral channels (right). 



OVERVIEW OF NULL RATE ANALYSIS 



Contribution 


Specification 


Best null (100 s) 


Null rate (7 h) 




10-* 


8.9 X 10-'5 


3.37 X 10-5 


Ns 


3.5 X 10-5 


1.7 X 10-6 


2.8 X 10-6 


Nx 


3.5 X 10-5 


4.6 X 10-"^ 


5.4 X 10-6 


Ne 


2 X 10-5 


3.9 X 10-6 


3.3 X 10-6 


NpoC 


10-5 


2.8 X 10-6 


2.2 X 10-5 


(JN,T=1s 


1.5 X 10-5 


3 X lO-'^ 


2 X 10-^ 


'^N,S,T=ls 


1.5 X 10-5 


3 X lO-'^ 


2 X 10-^ 


0'Af,e,T=ls 


3 X 10-6 


1.4 X 10-12 


1.5 X 10-" 


O'JV,T=100s 


10-5 




9 X 10-^ 



"estimation based on the other contributions 
Table 4. Overview of the analysis of the best null rate over 100 s and the null rate measured over 7 h, in average value 
and stability, and comparison with the original specifications of the bench. 



specifications of the bench. As we also demonstrate the capability of correcting injected disturbances, we should 
be able to achieve that result even with typical disturbances on the bench. 

6. CONCLUSION 

PERSEE showed a very high quality and performance well beyond the initial specifications. The large majority 
of goals has already been demonstrated. We have already presented some important lessons in Ref. |T31 about 
thermal stability and control loop feasibility. However, we can complete these remarks by other conclusions. 

Firstly, we elaborated a calibration scheme that can be applied to a space mission, although it must take 
into account the linear-parabolic drift. That procedure is fully automatic, and allows also to make an accurate 
analysis of the instrument. 

Secondly, we showed that a LQG controller can manage the disturbances of the OPD, to achieve a very deep 
null rate. An experimental sub-nanometer control has been demonstrated, even in the case of strong disturbances 
(more than 100 nm rms reduced to 1 nm rms). That result was very important for PERSEE, because we work in 
the near infrared, which is more sensitive to OPD variations than thermal infrared used in projects like DARWIN. 

Finally, we showed that thanks to a a very efficient design, we obtained null rates more than 10 times better 
than the initial conservative goal of lO-'', thus close to the requirement of a nuller for the observation of telluric 
planets. We achieved the deepest null ever achieved by a near-infrared nulling interferometer: 8.9 x 10-6 on a 
37% bandwidth, between 1.65 and 2.45 |^m. PERSEE bench shows also a very good stability over 100 s, and 
also over 7 hours. We still have to confirm this deep null rate over 10 hours, with the same type of stability, and 
with injected disturbances. 

The last steps will be the correction of linear-parabolic drifts, together with the use of a new source module, 
which includes a planet and an exozodiacal disk, presented in Ref. 15^ during this conference. 
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